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^ \ We derive the spin-dependent Bethe-Salpeter equation for qq bound state at finite temperature 



in the ladder and cross box diagram approximations. The imaginary time formalism is adopted in 
the adiabatic approximation where the flux tube thermal excitation is separated from the thermal 



excitation of the constituent quarks. The thermal excitation of the flux tube is implemented by 
Q \ using a temperature dependent potential inferred from the lattice gauge calculations. The thermal 



excitations of the constituent quarks are evaluated by summing Matsubara frequencies. The re- 



^ . sultant equation is used to study the thermal charmonium spontaneous dissociation mechanism at 

H ; 

5r , temperatures below the critical one of the phase transition to the quark gluon plasma. Our results 

show that when the thermal excitation of the constituent quarks is considered beside the flux tube 
thermal excitation, the spontaneous dissociation unlikely takes place in the hadronic phase. In 
contrary, when the thermal excitation of the flux tube is considered only in the calculations, the 
spontaneous dissociation likely takes place just below the critical temperature. 
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I. INTRODUCTION 



Charmonium suppression/enhancement in the heavy ions colhsion experiments has re- 
ceived much attention in the search of the quark gluon plasma (QGP) |l|, 0, |, ^, ^ j^. 
Whatever signal is chosen for the identification of QGP, contributions to that signal from 
conventional hadronic processes must be identified as backgrounds and removed from the 
data. Among of these signals is the signal from the charmoinum dissociation at temperatures 
below the QGP phase transition temperature Tc. There is a considerable uncertainty on the 
origin of the anomalous suppression due to the lack of a reliable information on J/ and xj 
dissociation below T^. It is known that J/ip production in hadron-hadron collisions is to a 
considerable extent due to the production and subsequent decay of higher excited cc states. 
Since different quarkonium states have different sizes, one expects that higher excited states 
will dissolve at smaller temperatures than the more tightly bound ground states^, ^. These 
facts may lead to a sequential suppression pattern in J/ ip yield in nucleus- nucleus collisions 
as a function of the energy density [||, |^. The charmonium can dissociate spontaneously 
when its mass spectrum becomes unbound and equal to a pair of final open charm mesons 
at temperatures below the critical temperature Tc. 

The spontaneous dissociation has been studied by using a temperature-dependent po- 
tential inferred from the lattice gauge calculations 0. In consideration of the spontaneous 
dissociation below T^, it is necessary to find the selection rules for the dissociation of a 
heavy quarkonium state with initial quantum numbers J, Li and Si into a pair of open 
charm mesons with a total spin S and a relative orbital angular momentum L. The total 
spin S is conserved. Parity conservation requires A = |L — Lj| = 1 in this dissociation. 
Hence J/ip and ip' dissociates into a pair of open charm mesons with L = 1, while x dissoci- 
ates into a pair of open charm mesons with L = or 2. These spin and angular momentum 
quantum numbers give rise to the selection rules for final meson states and alter the dis- 
sociation threshold energies and dissociation temperatures The spontaneous thermal 
dissociation also takes place when the charmonium absorbs a light meson such as a pion and 
their total mass becomes equal to a pair of open charm mesons via the threshold production 
interaction 

mj/^ + m.^ = mjj-^. (1) 
The selection rules still play a significant rule in this interaction. However, in the hadronic 
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phase the pion gas is accumulated in the media and it is hkely that charmonium dissociates 
thermally spontaneously by absorbing a pion. 

The other possibility is that the quarkonium can dissociate by colliding with light 
hadrons[|^, 0, ^, |TD[. This co-mover dissociation process must be understood and in- 
corporated in the simulation of heavy ion collisions before the QGP formation can be estab- 
lished through this signature. The nonrelativistic quark-interchange model has been used 
to evaluate the low-energy cross sections of heavy quarkonium in collision with light mesons 
such as TT, p and K in terms of wave functions and interactions at the quark level. It is 
found that as the temperature increases, the threshold energy decreases and the dissocia- 
tion cross section increases|l|, §]. Furthermore, any quantitative analysis of J /ip dissociation 
in nucleus-nucleus collision should include the effects of the meson mass modification in the 
dense and hot matter. 

The relativistic quasi-potential at finite temperature for the qq is essential to study the 
charm spectroscopy at finite temperature. Furthermore, it is expected that the relativistic 
corrections for the wave function shall be significant and they may change the results of the 
co-mover dissociation processes. The inter-quark temperature dependent potential modifies 
the spectrum significantly. Altogether, the thermal kinetic excitations of the constituent 
quarks are presumably supposed to change the full BS equation results. 

The outline of the paper is as follows. In Sec. II, we shall derive the BS equation at 
finite temperature. The adiabatic approximation is considered where the flux tube thermal 
excitation is separated from the kinetic thermal excitation of the constituent quarks. The 
interaction of the constituent quarks with the mean fields is also considered. In Sec. Ill, we 
study the thermal spontaneous dissociation at T < and give our conclusions. 

II. BETHE-SALPETER IN HOT AND DENSE MEDIUM 

The Bethe-Salpeter equation in the ladder approximation for the quark (1) and anti-quark 

(2) reads 

G^\P,)^{P)G',\P2) = ^ J 0^,nP.p'Mp')- (2) 

The quark propagator can be calculated by solving the Schwinger- Dyson equation in the 
ladder approximation. However, in the heavy-heavy and heavy-light meson bound states. 
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the self-energy corrections are small. Therefore, we consider the free quark propagator in 
the present calculations. When the meson is probed in the nuclear matter, the constituent 
quark interacts with the medium by coupling the scalar and vector mean fields. In the mean 
field approximations, the bound state equation probed in the nuclear matter becomes 

In symmetric isotropic nuclear matter, the vector field is simplified to time-like vector 
field. 



The bound state equation is simplified by introducing the following parameters. 

Ml = riiM* - 

M2 = 772M* - 

M12 = M* - [V^^^ + V^% 



(4) 



(5) 



The interaction with the vector mean fields can be absorbed in this new parameterization. 
Hence, the vector mean fields, effectively, couple with the bound state mass. This simphfies 
the equation drastically since the coupling with the vector fields will not appear any more. 
On the other hand, the scalar mean fields couple with the constituent quarks and modify 
their masses as follows. 



ml = mi — S^^^ 
The BS equation probed in the nuclear matter reads. 



(Mi+po)/3-p-7 + K 



^l^,^iP,p')W) 



where the effective constituent quark energy reads 



(po-M2)/3-p-7 + m; 
[ef{p)-{M2-poy-iS 



<(P) 



p2 + mf. 



(6) 



, (7) 



(8) 



The quasipotential equation in the nuclear matter has the same criteria of the normal quark- 
meson coupling model. Hence, the interaction with the mean scalar fields modifies the quark 



effective masses while the interaction with the vector mean fields modifies the effective bound 
state mass M*. 

When the bound state equation is embedded in the thermal bath, the constituent quarks 
are excited thermally as well as the flux tube. To simplify the equation, the thermal exci- 
tation of the flux tube between the constituent quarks is separated from the quark kinetic 
thermal excitations. In the adiabatic approximation, the flux tube is excited with the tem- 
perature independently of the thermal excitations of constituent quarks. Therefore, the flux 
tube excitation is not modifled with the thermal excitation of the quark propagator. We 
assume the instantaneous interaction approximation where the potential doesn't depend on 
the time-like component. Therefore, the kernel is approximated as follows, 

/C(p,p')^V(p-p',T,/x,). (9) 

The four-dimensional kernel integration is reduced to a three-dimensional one in the term 
of the instantaneous wave function 0(p) 

' I = ' J ^^^P - P'^T,^^,)<j>{p% (10) 

where the bound state wave function in the adiabatic and instantaneous approximations 
reads 

dpo. 



(p) = / ^m- (11) 

J ZTT 



In order to simplify the equation, we deflne the convention 

V(T,/i,)0 = I |^^V(p-p',T,/i,)</)(pO. (12) 

The three-dimensional bound state equation with the ladder graphs reads, 

^ . f dpo (Ml + po)f3 - 7 ■ p + m* (pp - Ma)/? - 7 ■ p + m* 

^^P^ = ' J ^(M,+,o)^-.f(p)+./(^' ^^^'^Wip)^^^^;^^^^ ^''^ 

The contribution of the cross box diagram justified by using the Eikonal approximation is 
added in the equation |TT|, |T^. It is given by 

J _ . f dpo (Ml + po)f3 - 7 ■ p + m* (po + M2)/3-7-p + m; 

^Cross - V 2vr (Ml + poY - ef (p) + ^6 ' ^''"^ efip) - (po + M^f - tS' ^ ^ 

This approximation is useful since it involves the same argument for both the ladder and 
cross-box graphs. It is noted that this term doesn't alter the invariance. This equation is 



symmetric in the particle labels and it reduces to the one body Dirac equation when either 
particle's mass becomes infinite. 

For the bound state problem in the center of mass frame M12 = M1 + M2 is the eigenvalue 
while the subsidiary condition Mi — M2 = {ml — m'^)/Mi2 should be used to define the 
energy difference. This condition guarantees that the correct nonrelativistic limit is used. 
The ladder and cross box diagrams are decomposed as follows, 



r(0) 



'Ladder -'Ladder ^ "^Ladder 



(1) 



and 



r ^ r(0) Al) 

^ Cross ^ Cross ^ Cross- 



(15) 



(16) 



The BS equation is simplified drastically when it is projected into the plane wave compo- 
nents. The positive and negative energy projectors read 

e:(p)±/3(7-p + m:) 



A2(P) = 
For the Ladder graphs we have 

r(0) 



Ladder Ladder Ladder 



2e|(p) 



(17) 



(18) 



where 



Ladder 



J 2n 



Af^(p) 



{po + Mi-et(p)+i5) 



(po - M2 + e|(p) - i(5) 



(19) 



^Ladder " ^ 



dpo 

271 



{po + Ml + el{p) - iS) 



Af (-P) 



(po - M2 - e*2{p) + iS) 



.(20) 



These terms lead to the Salpeter equation at zero temperature. The additional components 
in the ladder graphs read 



where 



Ladder 



. f dpo 



Ladder 



Af^(p) 



"^Ladder ~'~ "^Ladder' 



(21) 



Af^p) 



{po + Ml - €t{p) + iS) 



{po - M2 - €*2{p) + iS) 



(22) 
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j(— ) 

Ladder 



dpo 
27r 



A- (p) 



(23) 



{po + M, + el{p)-i5f ' '"^^""(po-M2 + e^(p)-i(5) 

These components vanish in the zero temperature case. The " + +" and " " which are 

essential to obtaining the Dirac hmit are absent at zero temperature unless the cross-box 
gluon exchange term is included. These terms correspond to the Z-graph contributions of 
time ordered perturbation theory. The contribution of the cross graphs read 



where 



^(++) 
Cross 



ji—) 
Cross 



dpo 
27r 

2tt 



Cross 



Af^(p) 



{po + Mi-el{p) + iS) 
P) 



J++ , J — 
Cross Cross' 



(po + Ml + €l{p) - i6) 



{po + M2 + e|(p) - i5) 
{po + M2 - e*2{p) + iS) 



The additional terms read 



Cross 



j(-+) 
Cross 



dpo 
2tt 



^(1) 
Cross 



Af^(p) 



{po + Mi-el{p) + i5) 



Cross Cross' 



A*}'\-p) 



dpo 

277 



AI^^^(P) 



-/3V(r,/.,)0(p)/3 



{po + M2-e*2{p) + i5) 



ipo + M, + elip)-t6r ' '"^"'"^"^"(po + A^2 + e^(p)-^^) 
These terms vanish at zero temperature by direct integration over residues. 
The bound state equation reads 



(24) 

(25) 
(26) 

(27) 

(28) 

(29) 



Hp) 



j(0) 

Ladder 



Cross box 



6Vr 



(1) 



T ' 



where 



j(i) 

Ladder 



j(i) 
Cross' 



(30) 



(31) 



Using the imaginary time formalism preforms the thermal extension. The right hand side is 
expanded with respect to the Matsubara frequencies 



271" I- Ladder Cross 



(1) 



00 

= iT E 



(0) 



J, 



(0) 



Ladder Cross 



6Vr 



+ 5V^^ 

P0=i0Jn -L 



= — i-^Res 
2 ^ 



(0) 



"Ladder Cross 



tanh(po/2r) + 



po 



(32) 



where a;„ = 7r(2n + 1)7. 

The quasipotential equation at finite temperature reads 

A«*(p)/3V(r,/x,)0/3AL')*(-p) 



Mi2-et(p)-e^(p) 

AL^^*(p)/3V(T,/x,)0/3Af*(-p) 
-Mi2-el(p)-e^(p) 



^^ + et(p)+e^(p) 



Mi2 

We have used the constraint (M^ — M|) = (mj^ — m^) in the cross-box diagram. The 
thermal distribution functions read 

U+-{T) = [1 - n(6t(p) - ml) - n(6;(p) - m^)] (34) 

N-+{T) = [1 - n(et(p) + ml) - n(e;(p) + m^)] (35) 

U++{T) = [1 - n(el(p) - ml) - n(e;(p) + m^)] (36) 

M-{T) = [1 - n(eKp) + ml) - n(e;(p) - m^)], (37) 

where n{x) — l/{exp{x/T) + 1) is the Fermi-Dirac distribution function. 
The approximation 

e*{p)±Mi^e*{p)±ml, (38) 

is justified by the confinement condition. In the absent of the mean fields, the confinement 

is satisfied by the absent of the time-like pole mi{p) + m2{p) > M12 or with a less constraint 



(33) 



condition \Jp'^ + mi{p) + ^Jp"^ + m^p) > M12 where the constituent quark mass mi{p) is 
obtained from the gap equation. Since, we don't solve the gap equation consistently with 
the BS equation but instead we use the free quark masses, this replacement is a reasonable 
one to avoid the time-like pole. 

The last term 5V!p^ that appears in the quasipotential equation is given by 

"^^^^(y) A(!)*(p)/jV(T,^,)<^/?Af*(-p) ^ ^^^A«*(p)/5V(T,/x,)<^/?aL^>(-p) 



5vf = 



+ 



Mi2 + (eUp)-e2(p)) ' ' Mi2+(eKp)-e^(p)) 

■ Ar(p)/^V(T,,,)0ML-)Vp)^^ AL-^p)/^V(T,..)0Mg)V^j 

. - (el(p) - e|(p)) - (^!(P) - 



where 



A^++(T) = [n{el{p)-ml)-n{e;{p) + m;)], 

M^4T) = [nielip) + ml) - n{e;{p) - m;)], 

M+-iT) = [n{el{p)-ml)-n{e;{p)-m;)], 

M^+iT) = [n(et(p)+mn-n(e;(p)+m;)]. 



(40) 
(41) 
(42) 
(43) 



The contribution of 6V^ vanishes at zero temperature and we have dropped it from the 
numerical calculations. 

The wavefunction in the plane wave components reads 

/ 



e*(p)+m* 



(44) 



V 4{p)+ml '/If; e'{p)+m* ej(p)+m* 



(P) 



g P 



V 



e*(p)+m* 

C(P) 



C(P) 



e*(p)+m* 

-C(P) 



C(P)^ 



e*(p)+m* 



e*(p)+m* 



P 



p 

V ei(p)+ml; 



-xip)7n^- 



^X(P) 



e*(p)+m* 

e*(p)+m* 



^ vfol ^ 



and 



c^P 



r^(p)z 



g P 



o^P 



e* (p)+m* Vf; e* (p)+m* e* (p)+m* 

^(p)z(S^ ^(p 



T^(p) 



(45) 



(46) 



(47) 



Hence the BS equation is decomposed into the four independent radial plane waves, 



^ v{p) ^ 






-(T)v+r 


A/; 




' Ar+ 




A/; 


_(T)V+^^ 






ix{p) 




A/'+ 




A/'+ 




^ Ar+ 




A/'+ 


4T)V++ 




^x(q) 


#(p) 


AT, 


4T)vir 




-(T)V± 


' AT- 


-{T)VZZ 




4T)vr± 




i^(q) 


I CP) J 




U- 


4T)V1+ 






' M- 










I Cq) j 



(48) 



where the effective propagator reads 



Q = 











-1 

e'+e;-(mf -m*^)/Mi2 











1 

eJ+€*+(m*^-m*^)/Mi2 










e*+e|+Mi2 y 



(49) 
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In the present calculations, we have adopted the kernel interaction as the vector Coulomb 
and scalar linear confining potentials 



V(T,/x,) = 7;. X rVy{T) + 1 x lyg(T) 



(50) 



where Vy and Vg represents the short and long range interaction, respectively. The tem- 
perature dependent potential for the long-range interaction has been computed H. It reads 



Sn — 



TT 

12 



6 



arctan(2rT) 



+ 



ai{T) 



TT 2 2 

— T H — T arctan , 
3 3 V2rT 



r H m 

2 



1 + {2rTy 



(51) 



The string tension cr(T) includes thermal corrections to the zero-temperature string tension. 
It is assumed to have the form 



/ ^2X1/2 



(52) 



for temperatures below T^, the critical temperature of the deconfinement phase. It is urged 
that the string tension does not vanish at the critical temperature [Q. In the numerical 
calculations, we have taken ctq = 0.22GeV^ and b = 0.99. The short-range interaction is 
taken as 



a 



Vy{T) = exp(-/iscreen(T)r). 



(53) 



Since the Debye screening becomes active only in the deconfined phase, we set /iscreen = 
for T < Tc. The fitting parameters [a„ So] are taken [0.23, 140MeV] and [0.27, -156MeV] for 
charmonium and a pair of open charm mesons, respectively. The constituent quark masses 
are taken as 300 and 1300 MeV for and rriu^d and mc, respectively. 



III. DISCUSSION AND CONCLUSIONS 



We have derived the relativistic bound state equation for qq probed in a hot and dense 
matter. It is found that the scalar mean fields modify the constituent quark mass by m* = 

~ J2s QqsO's while the vector mean fields modify the effective bound state mass by M*^ = 
Mgg + J2v Qmv^v where Mqq is the bound state mass without vector main fields. Therefore, 
the vector mean fields decouple from the constituent quarks in the relativistic equation while 
the scalar mean fields modify the relativistic bound state equation. 
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The adiabatic approximation is considered in the derivation of the BS equation probed 
in the hot bath. The thermal excitations of the constituent quarks are separated from the 
thermal excitations of the flux tube. The flux tube thermal excitation is taken into account 
by adopting the central inter-quark potential derived by the lattice gauge calculations [Q. 
However, the constituent quark thermal excitation is considered in the calculation by deriv- 
ing the BS equation in the ladder and cross-box approximations using the imaginary time 
formalism where the summation of the Matsubara frequencies introduced in the constituent 
quark propagator is evaluated. Since we have considered the free quark propagator instead 
of the exact gap equation in BS equation, we adopt the approximation given by Eq. ^ 
in the Fermi-Dirac distribution functions to justify the confinement condition. This condi- 
tion ensures the confinement and avoids the thermal ionization of the bound state. We have 
only considered the terms of the ladder and cross-box diagrams which have zero temperature 



counterparts |jT3|. The other terms that appear in the ladder and cross-box diagrams at finite 
temperature but vanished at zero temperature are dropped from the calculations. However, 
the contribution of these terms is negligible. We are not sure if these terms are cancelled by 
counterparts that come from other diagrams with a thermal origin. The BS equation with 
the quantum numbers {JLS} is solved exactly. The full relativistic quark kinematics and 
their thermal excitations are considered exactly in the calculations. These fully relativistic 
kinetic effects become significant at the energies involved in relativistic heavy ions collisions. 

We have solved exactly the BS equation where both the constituent quarks and flux tube 
thermal excitations are considered in the calculations. The charmonia and charmed meson 
pair masses versus temperature are displayed in Fig. 1. The meson masses decrease rapidly 
at first with temperature but when the temperature reaches about 50 MeV this decrease is 
slowed down considerably. The mass differences between the various charmonium spectra 
with different quantum numbers decrease drastically with the temperature and become very 
small at high temperatures. The hyperfine splitting for the D-meson decreases significantly 
and becomes negligible when the temperature exceeds T ~ 100 MeV. Near but below the 
critical temperature the threshold production energy difference between the pair of open 
charm mesons and charmonium is roughly of the order 3-4 pions. The mass difference 
between Xc and a pair of open charm D{^Si) is about 612 MeV while the energy difference 
between the J/ip (or 77^) and a pair of open charm D{^Si) (or D{^Sq)) is almost 613 MeV. 
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Just below the critical temperature, we still have 

Mj,-^ » M^^ + M^^ Mj/^ + M^. (54) 

Therefore, although the J/ip can be easily excited thermally to higher angular momentum 
states, it and its excited states still can't reach the threshold energy production of the open 
charm mesons pair even if it absorbs a pion from the media. In the heavy ions collisions 
charmonium is surrounded by the pion gas. 

We have also considered the BS equation with only the flux tube thermal excitations. 
In this approximation the constituent quarks are assumed heavy and cannot be thermally 
excited. The bound state equation simply reduces to the zero temperature BS equation but 
with a temperature dependent potential. In this limit, the quark and anti-quark partition 
functions reduce to A/±± — A4±± — 1. The meson spectra masses versus temperature are 
displayed in Fig. 2. The charmonia and the open charm pair masses decrease slowly with 
the temperature. However, their masses start to drop significantly just below the critical 
temperature. Furthermore, the energy differences between states with different quantum 
numbers does not change much with the temperature. Therefore, the selection rules persist 
until almost the critical temperature. The mass difference between Xc and a pair of open 
charm D{^Si) is found to be 52, 100.1 and 192 MeV for Xo,Xi and X2, respectively. The 
masses difference between J/t/j and a pair of open charm D{^Si) is found to be 147 MeV 
while it is 158 MeV for rjc and a pair of open charm D(^So)- The minimum energy difference 
between the pair of open charm mesons and charmonium becomes of the order of a pion 
mass or less below the critical temperature. The threshold production energy is given as 

M^^ < M^, + M^^ Mj,^ + M^. (55) 

Moreover, the production energy is of the same order as the mass fitting deviation for Xc- 

^xo - ^DD ^ W^f^ - Mfit|. (56) 

Therefore, the J/ip can be thermally excited to a higher orbital angular momentum state 
and then decays to a pair of open charm mesons either spontaneously or by absorbing a 
pion and then decays to open charm mesons pair spontaneously. This interaction raises 
drastically the J/%1) or Xc suppression probability. We have compared the ground state 
masses for charmonia and open charm mesons pair with and without constituent quark 
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kinetic thermal excitations. The results are displayed in Fig. 3. The hyperfine splitting 
becomes less important when the constituent quark thermal excitations are involved in the 
calculations and the threshold production energy to a pair of open charm mesons becomes 
much higher. 

In the conclusion, our results show that there is no thermal route for the spontaneous 
thermal J/ip suppression in the hadronic phase when the thermal excitation of the con- 
stituent quarks is taken into account beside the temperature dependent potential derived 
from the lattice gauge calculations. 
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FIG. 1: Charmonium and a pair of open charm mesons DD spectra versus temperature with 
thermal flux tube and quark excitations. 
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FIG. 2: Charmonium and a pair of open charm mesons DD spectra versus temperature with only 
thermal flux tube excitation. 



15 



4.4 
4.2 
4.0 
3.8 

> 

O 3.6 

o 

(0 3.4 

</> 

CO 

2 3.2 

3.0 
2.8 
2.6 



— % (I) 
J/*(i) 

D(^So) pair (I) 

D('Si) pair (I) 

% (II) 

- — J/^(ll) 

D(^So) pair (II) 

- - D('Si) pair (II) 




0.0 



0.05 



0.1 



0.15 



0.2 



T( GeV) 



0.25 



FIG. 3: Comparison between the spectra of charmonium and a pair of open charm mesons pro- 
duction with the thermal flux tube excitation and the thermal quark and flux tube excitations. 
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